RANDOM WALKS OF CIRCLE PACKINGS 



TOMASZ DUBEJKO 

Abstract. A notion of random walks for circle packings is introduced. The geom- 
etry behind this notion is discussed, together with some applications. In particular, 
we obtain a short proof of a result regarding the type problem for circle packings, 
which shows that the type of a circle packing is closely related to the type of its 
tangency graph. 



1. Introduction 

Circle packings have been introduced by Thurston to build finite Riemann map- 
pings and to conjecture the Finite Riemann Mapping Theorem ([Tl], [T2]). They 
have been used to construct discrete analytic function theory ([BStl], [Dl], [D2], 
[HR], [RS], [Stl]). Interesting connections between circle packings and graph- 
theoretical aspects of their tangency patterns have been investigated in [Stl], and 
recently in [Mc], [BeSc], [D3], and [HSc2]. 

Here we are interested in developing geometric notion of random walks and 
harmonic functions for circle packings. This notion will be first described for circle 
packings with mutually disjoint interiors, i.e. univalent packings, and then will 
be extended to branched and overlapping-angles circle packings. Our construction 
will yield explicit harmonic functions for random walks of packings, which are also 
Dirichlet-finite when associated packings have compact carriers. It will also follow 
that the real and imaginary parts of finite Riemann mappings of [Tl] are harmonic 
for underlying random walks. Moreover, we obtain short (and elementary, as they 
can be) proofs of the following. 

Theorem 1.1. Let IK be a triangulation of an open disc in the complex plain C. 
If the simple random walk on the 1-skeleton K 1 of K is recurrent then there exists 
a univalent circle packing for K whose carrier is C. Conversely, z/K is of bounded 
degree and there exists a univalent circle packing for K whose carrier is the unit 
disc D then the simple random walk on K 1 is transient. 

Theorem 1.2. LetK. be a triangulation of an open disc. IfK is of bounded degree 
and transient then there exist non-constant, Dirichlet-finite, harmonic functions for 
the simple random walk on IK 1 . 

Several remarks regarding the above theorems should be made. It was shown 
([BStl], [HScl]) that given K as in Theorem 1.1 there exists a univalent circle 
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packing for K whose carrier is C or D, but not both. The problem of deciding 
whether the carrier is C or D has become to be known as the type problem for 
circle packings (see [BStl]), by analogy to the classical type problem for graphs. 
The first answer regarding this problem was given in [Mc] for IK of bounded degree. 
Theorem 1.1 subsume the result of [Mc], and has been very recently established in 
[HSc2] using interesting but involved combinatorial extremal length techniques. We 
give here simple and independent proof of this result as an example of an application 
of our new notions. 

The advantage of working with random walks of packings, rather then with sim- 
ple random walks of their tangency graphs, is that there exist non-trivial harmonic 
functions associated with them which have a geometric interpretation. Using geo- 
metric properties of these functions and the fact [Me] that the existence of Dirichlet- 
finite harmonic functions is preserved by rough isometries of networks, Theorem 1.2 
easily follows from Theorem 1.1. We note that Theorem 1.2 was originally proved 
in [De] using cohomology theory and also follows from a new result of [BeSc]. We 
also observe that Theorem 1.2 and [So2]/[Me] imply that there exist non-constant, 
Dirichlet-finite, harmonic functions for graphs of bounded degree which are roughly 
isometric to the 1-skeleton of a triangulation of Theorem 1.2; this includes, for ex- 
ample, planar trees of bounded degree and planar cell decompositions of bounded 
degree with a uniform bound on the number of edges in each cell. 

Acknowledgments. The author would like to thank K. Stephenson for interesting 
discussions and P. M. Soardi for the reference [Me]. Figure 2 was created with the 
help of Stephenson's software CirclePack. 

2. Preliminaries 

Since there are number of publications on the subject of circle packings, we 
skip formal definitions. Our terminology will be consistent with that of [Dl] and 
[DSt] (see also [BStl]) for circle packings without overlapping-angles. The reader 
is refered to [BoSt] and [T2] for information about circle packings with angles of 
overlaps. 

We will now set up our notation. If IK is a triangulation, finite or infinite, with or 
without boundary, of a disc in C then we will write K° and IK 1 for sets of vertices 
and edges of IK, respectively. If V is a circle packing for IK then C-p{v) will denote 
the circle in V corresponding to a vertex v G K°, r-p(v) will be the radius of C-p(v), 
and S-p : K — > carr("P) will be the simplicial isomorphism of K with the carrier 
carr(P) of V determined by: S-p(v) := the center of C-p(v), v G 1K°. 

As far as the subject of random walks is concerned, the survey article [W] is 
more than a sufficient source of information for our purposes (see also [Sol]). 

If K is as above then (K,p) will denote a random walk on the 1-skeleton of K 
(i.e., on the graph IK 1 ), where p : K° x K° — > [0, 1] is the transition probability 
function. To simplify the notation, if tt, v G are adjacent then we will write 
u ~ v, and we will denote their common edge by uv. In the sequence we will be 
working with random walks given by conductances. If c : K° x K° — > (0, oo) is such 
that 

(1) c(u, v) > iff u ~ v, 

(2) c(u,v) = c(v,u), 

(o\ „/„. „A n £ — „n „. wO 



RANDOM WALKS OF CIRCLE PACKINGS 



3 



then the random walk (K,p) given by p(n, v) := c(u,v)/(^2 wr ^ u c(u,w)) will be 
called revesible with the edge conductance between u and v, u ~ v, equal to c(u, v). 
The simple random walk on IK 1 is a random walk with all edge conductances equal 
to 1. A real- valued function / : K° — > R is said to be harmonic for (IK,p) if 

f(v) = u)f(u) for every v G intIK , 

where int K is the set of interior vertices of IK. A function / : K° — > R is said to be 
Dirichlet-finite for a reversible random walk (K, p) with edge conductances c(-u, u) 
if 

(/(«)-/(«)) c(u,v)<oo. 

A random walk (K, p) is called recurrent if probability of visiting any vertex infin- 
itely many times is 1; otherwise is called transient. Finally, a graph K 1 is called 
recurrent {transient) if the simple random walk on IK 1 is recurrent (transient). 

Two important results linking the existence of harmonic functions with the type 
of a random walk are summarized in the following. 

Theorem 2.1. [W, Thm. 4.1, Cor. 4.15] Let K be a triangulation of an open disc 
in C. 

(1) If(K,p) is recurrent then there are no non-constant, non-negative harmonic 
functions for (IK, p). 

(2) If (IK, pi) and (IK, p 2 ) are two reversible random walks with edge conduc- 
tences c\ and c^, respectively, such that c\{u,v) > C2(u,v) for all uv G IK 1 
then the recurrence o/ (IK, pi) implies the recurrence o/(lK, p 2 ). 

We also need to recall the notion of rough isometry. Two metric spaces (X, dx) 
and (Y, dy) are said to be roughly isometric if there exist a map ip : X — > Y and 
constants A > and B > such that 

A^dxix, x')-B < d Y (ip(x), ip(x')) < Ad x (x, x') + B for all x, x G X, 

and 

d Y (y,ip(X)) < B for every y G Y. 

For a graph X and a reversible random walk (X, p) on X with edge conductances 
c(-, ■), we define the associated metric space (X, d p ) as follows. First we define the 
length of an edge path 7 in X by 

xx'e-y v ' 

Then the metric d v on X is given by 

d p (x,x') = inf{£(7) : 7 is a path in X with end points x and x'}. 

1 

We will say that two reversible random walks (X, p) and (Y, q) are roughly iso- 
metric if (X, d p ) and (Y, d q ) are roughly isometric. We finish preliminaries with the 

-f„i !„„.:„„. „..1j- -tu„-t u„„ 1 „u„„.„ ri\/r„l /„ — „i„„ re* ol\ 
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Theorem 2.2. [Me, Thm. 3.6] Let X and Y be two graphs of bounded degree. 
Suppose (X, p) and (Y, q) are roughly isometric reversible random walks such that 
Cp(x,x') < k and c q (y,y') < k for some constant k and all edges xx' and yy' , 
where c p and c q are edge conductances of p and q, respectively. Then there exists a 
non-constant, Dirichlet- finite, hermonic function for (X, p) iff there exists one for 
(Y,g). 

3. Random walks of univalent packings 

In this section we will be working only with univalent circle packings, i.e. circle 
packings whose circles have mutually disjoint interiors. Suppose V is a univalent 
packing for IK, and v G int K° . If u is a neighbor of v then there are exactly two ver- 
tices w' and w" that are adjacent to both u and v. Let r!p(uv) and r'^{uv) be radii of 
two orthogonal circles to triples C-p(u) ,C-p(v) ,C-p(w') and C-p(u),C-p(v),C-p(w"), 
respectively (Fig. 1). 




Figure 1. Orthogonal circles. 



We define the edge conductance between u and v induced by V by 

/*\ i \ r' v (uv) + r'^{uv) 

r v (u)+r r (v) 

We can now introduce 

Definition 3.1. If V is a univalent circle packing for K then the random walk 
(K,pp) induced by V is defined as follows: 

f cr(u,v) — if v & intK , u ~ v, and u ^ v 
0, if v e int K° and u = v,orv £ int K° and u ^ v 
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Remark 3.2. 

1) If IK has no boundary then (K,pp) is reversible with edge conductances 
c-p(tt, v). If IK has boundary then the boundary is an absorbing barrier for 
{K,p v ). 

2) From Herron's formula it follows that c-p(u,v) can be written explicitly in 
terms of radii of circles of V, 

( \ _ ^/rv{u)r v (v) ( I rv (w') I r v (w") 

^V\u, U) r v (u)+r-p(v) I V r v (u)+r P (v)+r v (w') V r v (u)+r v (v)+r v (w") 



where w' and w" are vertices adjacent to both u and v. 
3) If IK has no boundary and is of bounded degree, and V is a univalent pack- 
ing for IK then the Ring Lemma [RS] implies that there is a constant k, 
depending only on the degree of IK, such that c-p(u,v) G for all 

edges uv G IK 1 . In particular, by Theorem 2.1.(2), (IK, p-p) is recurrent iff 
IK 1 is recurrent. Moreover, (K, p-p) and the simple random walk on IK 1 are 
roughly isometric. 

We will now show that with every random walk induced be a circle packing there 
are associated two non-trivial harmonic functions. 

Lemma 3.3. Suppose V is a univalent circle packing for K and (IK, p-p) is the 
induced random walk on IK 1 . Let 9?-p : IK — > R and Ss-p : IK — > R be defined 
by $l-p(v) := Re[S-p(v)) and $s-p(v) := Im^-p (?;)), respectively, where S-p : K — > 
carr(P) C C is the associated simplicial isomorphism, and Re and Im denote real 
and imaginary parts of a complex number, respectively. Then functions Sft-p and Q-p 
are harmonic for (IK, p-p). 

Before we give a proof of the lemma we need the following proposition. 

Proposition 3.4. Let zo, . . . ,z m G C and let 7 = zqZi U Z1Z2 U ■ • • U z m zo be a 
closed, piecewise-linear, oriented path in C built out of oriented segments ZjZj + i 
(j mod (m+ 1)). Denote by rjj the unit normal vector to ZjZj+i. Then 

m 

J2\ z j+i ~ z j\Vj = 0, 
3=0 

where \z\ is the modulus of z. 

Proof. Write z j+1 - Zj = \z j+1 - Zj\e iLp ^ where ipj G [0, 2n). Then rjj = e i ^ +7T/2 \ 
Thus 

mm m 

= e-/ 2 ($> J+1 - Zj )) = Y^e^Zj+i - zj\e^ = £ \z j+1 - Zj \ Vj . □ 

j=0 j=0 j=0 



Proof of Lemma 3. 3. It is sufficient to show that 

cp(v, u)(Sp(u) — Sj>(v)) = for every v G intK . 
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For v E int K° let v°, . . . , v m be its heighboring vertices listed in the positive orien- 
tation. Let Zj be the center of the orthogonal circle to the triple of circles C-p(v), 
C-p(v j ), and Cj>(v j+1 ). From (*) and Proposition 3.4 it follows 

^2 c r(v, u) (S v (u) - S v (v)) = ^2 Cr p( v i vJ ) (Sv(v J ) - S r (v)) 

E\ Z 3 - Z 3-l\ (c ( 3\ C I \\ I \Sv(v j ) - S-p(v) 
— (S-p{V J ) - O-plV)) = > \Zj - Zj-i\ — -— = 0. 
r v (v') +r v (v) v v ' K " ^ 1 3 3 r v (v3)+r P (v) 

□ 

Suppose that V and Q are two univalent circle packings for K. Then the circle 
packing map <SV,Q := Sq ° S-p 1 '• carr(7 : ') — > carr(Q) maps centers of circles in 
V to centers of the corresponding circles in Q, and can be regarded as a discrete 
analog of an analytic function (see [BStl], [Dl], [RS], [Tl]). In particular, functions 
3?q and ^sq can be viewed as a harmonic function and its conjugate for (K,pq), 
respectively. 

In [RS] (see also [Stl], [HR]) it was proved that finite Riemann mappings (which 
are circle packing maps) of [Tl] can be used to approximate the classical Riemann 
mapping. If O is a Jordan domain in C, A C O is compact, and e > then 
for every sufficiently small mesh regular hexagonal circle packing V filling in Q 
there exists combinatorially equivalent (i.e., with the same tangency graph) circle 
packing Q filling in the unit disc D such that \Sj>^q{z) — t(z)\ < e for all z e A, 
where r is the Riemann mapping from O onto D (Fig. 2). Using the fact that 
the ratio function Sip q : carr(P) — > R, Sip q(v) := vq(v) /r-p(v), associated with 
S-p,Q approximates \r'\ in A, it is not hard to show that the harmonic measure of 
dA is approximated by escape (hitting) probabilities of (K, Pq)\k a , where IK is the 
tangency graph of V (and Q), Ka is the maximal simply connected 2-subcomplex of 
IK such that {C-p(v)} veK o^ C A, and (&,Pq)\k a is the random walk on Ka obtained 
by restricting (K,pg) to IK^ with 8Ka as absorbing barrier. It is our conjecture 
that the harmonic measure of dVt is approximated by escape probabilities of (IK, pg). 
However, this problem seems to be rather difficult to settle because up-to-date 
results show that r and \r'\ can be approximated uniformly only on compacta of Vt 
by circle packing maps. 
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We will now give proofs of Theorems 1.1 and 1.2. 

Proof of Theorem 1.1. Suppose IK 1 is recurrent. From [HScl, Thm. 0.3] there exists 
either a univalent circle packing for K with carrier equal to C or a univalent circle 
packing with carrier equal to D, but not both. Denote this packing by V. From 
Remark 3.2.2) it follows that 



Since edge conductances of the simple random walk on IK 1 are all equal to 1, Theo- 
rem 2.1.(2) implies that (K, p-p) is recurrent. Now, if carr^) were equal to D then, 
by Lemma 3.3, there would exist a bounded, non-constant, harmonic function for 
(K,pp). This, however, is impossible by Theorem 2.1.(1). Hence carr('P) = C. 

Suppose now that K is of bounded degree and there exists a univalent circle 
packing V for K such that carr('P) = D. To show that IK 1 is transient it is sufficient, 
by Remark 3.2.3), to show that (K, p-p) is transient. This follows immediately from 
Lemma 3.3 and Theorem 2.1.(1). □ 

Remark 3.5. The result in Theorem 1.1 is best possible. It is easy to construct (see 
[BeSc]) a circle packing with carrier equal to C whose tangency graph is transient 
(and, of course, of unbounded degree). 

Proof of Theorem 1.2. From Theorem 1.1 and [HScl, Thm. 0.3] there exists a 
univalent circle packing V for IK with carr('P) = D. Let $lp be as in Lemma 3.3. 
Let M be the degree of K 1 . Then, by Remark 3.2.2), 



Hence, $tp is a non-constant, Dirichlet-finite, harmonic function for (IK, p-p). By 
Theorem 2.2 and Remark 3.2.3), it follows that IK has non-constant, Dirichlet-finite, 
harmonic functions. □ 

From the last proof we obtain the following corollary. 

Corollary 3.6. Suppose IK is a triangulation of an open disc in C. If IK 1 is of 
bounded degree and transient, there exists a circle packing V for IK with carr('P) 
compact. Moreover, functions 3?-p and Qp are non- constant, bounded, Dirichlet- 
finite, and harmonic for the random walk (IK, p-p). 



We will generalize the notion of random walks induced by circle packings to 
non-univalent circle packings. For branched packings without angles of overlaps 
(see [Dl], [D2], [DSt]) one can use Definition 3.1. However, in order to have a 
comprehensive approach, we proceed as follows. 

Suppose V is a packing for K, possibly branched and with prescribed angles of 
overlaps (see [BoSt], [T2]). Let v G intIK . As before, if u is a neighbor of v then 




(9Mu) -^ v (v)) 2 c v {u,v) < k 0M«)-3Mf)) 2 




<k ^ (rr{u) +r v (v)) 2 < k Y 2 ( r r( u ) + r r( v )) < 2k Y1 Mr r( v ) < 2kM - 




4. Random walks of non-univalent circle packings 
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Figure 3. Radical centers and the edge conductance. 

z' and z" be the radical centers (see [C], [Y]) of triples of circles Cp(-u), C-p(v), 
C-p{w') and C-p(w), Cp(i'), C-p(w"), respectively. 

We set the edge conductance between v and u induced by V to be 



c P (u,v) 



\z'-z"\ 



\S v (u) - S v (u)\ 
We are now ready to describe the random walk of V. 

Definition 4.1. If V is a circle packing for IK, possibly branched and with pre- 
scribed angles of overlaps, then the random walk (K, p-p) induced by V is defined 
as follows: 

{ cr(u,v) — •£ intK , u ~ v, and u ^ v 
0, if v G int K° and u = v,oxv i int K° and u ± v 

1, if v int K and u = v. 

Remark \.2. If K has no boundary then it follows from the definition that (K,pp) 
is reversible with edge conductances c-p{u,v). If all angles of overlaps in V are 
(i.e., V is a locally univalent or branched circle packing) then c-p(u, v) are given by 
the formula in Remark 3.2.2). 

We finish with the following generalization of Lemma 3.3; the proof is similar to 
that of Lemma 3.3 and is left to the reader. 

Lemma 4.3. Suppose V is a circle packing for K, possibly branched and with 
prescribed angles of overlaps. Then both functions 3?p, 3p : K° — > R, Kp(w) := 
Re(S-p(v)) , ^s-p( v ) '■= ^- m {Sv( v )) , are harmonic for (K, p-p). 
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